Minimal Polynomial
Jordan Chains
Let A be an nxn matrix and Ā be its Jordan canonical form.
A ˉ = B − 1 A B , where B is invertible and composed of the basis vectors for the N ( A − λ i I ) m i \bar A = B^{-1}AB \text{, where B is invertible and composed of the basis vectors for the } N(A-\lambda_iI)^{m_i} A ˉ = B − 1 A B , where B is invertible and composed of the basis vectors for the N ( A − λ i I ) m i
rank ( A ) = rank ( B A ) = rank ( A B ) = rank ( A ˉ ) \text{rank}(A) = \text{rank}(B A) = \text{rank}( A B) = \text{rank}(\bar A) rank ( A ) = rank ( B A ) = rank ( A B ) = rank ( A ˉ )
dim ( N ( A − λ i I ) ) = dim ( N ( A ˉ − λ i I ) ) \text{dim}(N(A-\lambda_iI)) = \text{dim}(N(\bar A - \lambda_iI)) dim ( N ( A − λ i I )) = dim ( N ( A ˉ − λ i I ))
A ∈ C n × n A ˉ = J = B − 1 A B C n = N ( A − λ 1 I ) m 1 ⊕ ⋯ ⊕ N ( A − λ k I ) m k We look for basis vectors for N ( ( A − λ i I ) m i ) B = [ B 1 B 2 ⋯ B n ] B i = [ b i 2 b i 3 ⋯ b i r i ] where b i , ⋯ , b i r i have to be basis vectors for N ( ( A − λ i I ) m i ) A \in \mathbb{C}^{n \times n} \\
\bar A = J = B^{-1}AB \\
\mathbb{C}^n = N(A-\lambda_1I)^{m_1} \oplus \cdots \oplus N(A-\lambda_kI)^{m_k} \\
\ \ \ \ \ \ \text{We look for basis vectors for } N((A - \lambda_iI)^{m_i}) \\
B = \begin{bmatrix} B_1 & B_2 & \cdots & B_n \end{bmatrix} \\
B_i = \begin{bmatrix} b_i^2 & b_i^3 & \cdots & b_i^{r_i} \end{bmatrix} \text{ where } b_i, \cdots, b_i^{r_i} \text{ have to be basis vectors for } N((A-\lambda_iI)^{m_i}) A ∈ C n × n A ˉ = J = B − 1 A B C n = N ( A − λ 1 I ) m 1 ⊕ ⋯ ⊕ N ( A − λ k I ) m k We look for basis vectors for N (( A − λ i I ) m i ) B = [ B 1 B 2 ⋯ B n ] B i = [ b i 2 b i 3 ⋯ b i r i ] where b i , ⋯ , b i r i have to be basis vectors for N (( A − λ i I ) m i )
Let M i : = A − λ i I M_i := A - \lambda_iI M i := A − λ i I , and let's choose a vector x such that x ∈ N ( M i j ) x \in N(M_i^j) x ∈ N ( M i j ) or M i j x = 0 M_i^jx = 0 M i j x = 0 , but x ∉ N ( M i j − 1 ) x \notin N(M_i^{j-1}) x ∈ / N ( M i j − 1 )
Now consider the chain of vectors x , M i x , M i 2 x , ⋯ , M i j − 1 x x, M_ix, M_i^2x, \cdots, M_i^{j-1}x x , M i x , M i 2 x , ⋯ , M i j − 1 x , which is linearly independent. We can see this by contradiction. Suppose that the chain is linearly dependent. Then there exists a vector M i k x M_i^kx M i k x such that M i k x = ∑ l = 0 k − 1 c l M i l x M_i^kx = \sum_{l=0}^{k-1} c_lM_i^lx M i k x = ∑ l = 0 k − 1 c l M i l x , where c l ∈ C c_l \in \mathbb{C} c l ∈ C . But then M i k x = 0 M_i^kx = 0 M i k x = 0 , which is a contradiction. Thus, the chain is linearly independent.
{ M i m i − 1 x , M i m i − 2 x , ⋯ , M i x , x } is linearly independent \{ M_i^{m_i-1}x, M_i^{m_i-2}x, \cdots, M_ix, x \} \text{ is linearly independent} { M i m i − 1 x , M i m i − 2 x , ⋯ , M i x , x } is linearly independent
┌──────────────────────────────────┬──────────┐
│ │ N(Σᵢᵐ) │
│ - x └──────────┤
│ ┌───────────────────────────┬──────────┤
│ │ │ N(Σᵢᵐ⁻¹) │
│ │ - Σᵢx └──────────┤
│ │ ┌────────────────────┬──────────┤
│ │ │ │ N(Σᵢ²) │
│ │ │ - Σᵢᴹᶦ⁻² └──────────┤
│ │ │ ┌───────────────┬──────────┤
│ │ │ │ │ N(Σᵢ) │
│ │ │ │ └──────────┤
│ │ │ │ - Σᵢᴹᶦ⁻¹ │
└──────┴──────┴────┴──────────────────────────┘
Example : Find the Jordan canonical form of A = [ 1 1 − 1 − 1 3 2 0 0 1 ] A = \begin{bmatrix} 1 & 1 & -1 \\ -1 & 3 & 2 \\ 0 & 0 & 1 \end{bmatrix} A = ⎣ ⎡ 1 − 1 0 1 3 0 − 1 2 1 ⎦ ⎤
Solution :
d ( s ) = det ( s I − A ) = ( s − 1 ) [ ( s − 3 ) ( s − 1 ) + 1 ] = ( s − 2 ) 2 ( s − 1 ) d(s) = \text{det}(sI - A) = (s-1)[(s-3)(s-1)+1] = (s-2)^2(s-1) d ( s ) = det ( s I − A ) = ( s − 1 ) [( s − 3 ) ( s − 1 ) + 1 ] = ( s − 2 ) 2 ( s − 1 )
λ 1 = 2 , λ 2 = 1 \lambda_1 = 2 \text{, } \lambda_2 = 1 λ 1 = 2 , λ 2 = 1
Σ 1 = A − 2 I = [ − 1 1 − 1 − 1 1 2 0 0 − 1 ] \Sigma_1 \ = \ \ A-2I \ \ = \begin{bmatrix} -1 & 1 & -1 \\ -1 & 1 & 2 \\ 0 & 0 & -1 \end{bmatrix} Σ 1 = A − 2 I = ⎣ ⎡ − 1 − 1 0 1 1 0 − 1 2 − 1 ⎦ ⎤
dim N ( Σ 1 ) = 3 − 2 = 1 ≠ r 1 = 2 \text{ dim }N ( \Sigma_1) = 3-2 = 1 \neq r_1 = 2 dim N ( Σ 1 ) = 3 − 2 = 1 = r 1 = 2
Σ 1 2 = ( A − 2 I ) 2 = [ − 1 1 − 1 − 1 1 2 0 0 − 1 ] [ − 1 1 − 1 − 1 1 2 0 0 − 1 ] = [ 0 0 4 0 0 1 0 0 1 ] dim N ( Σ 1 2 ) = 3 − 1 = 2 = r 1 \Sigma_1^2 = (A-2I)^2 = \begin{bmatrix} -1 & 1 & -1 \\ -1 & 1 & 2 \\ 0 & 0 & -1 \end{bmatrix} \begin{bmatrix} -1 & 1 & -1 \\ -1 & 1 & 2 \\ 0 & 0 & -1 \end{bmatrix} = \begin{bmatrix} 0 & 0 & 4 \\ 0 & 0 & 1 \\ 0 & 0 & 1 \end{bmatrix} \text{ dim }N(\Sigma_1^2) = 3-1 = 2 = r_1 Σ 1 2 = ( A − 2 I ) 2 = ⎣ ⎡ − 1 − 1 0 1 1 0 − 1 2 − 1 ⎦ ⎤ ⎣ ⎡ − 1 − 1 0 1 1 0 − 1 2 − 1 ⎦ ⎤ = ⎣ ⎡ 0 0 0 0 0 0 4 1 1 ⎦ ⎤ dim N ( Σ 1 2 ) = 3 − 1 = 2 = r 1
No need to check Σ 2 \Sigma_2 Σ 2 because r 2 = 1 r_2 = 1 r 2 = 1
m ( s ) = ( s − 2 ) 2 ( s − 1 ) = d ( s ) m(s) = (s-2)^2(s-1) = d(s) m ( s ) = ( s − 2 ) 2 ( s − 1 ) = d ( s )
From the minimal polynomial, we can see that the Jordan canonical form will have a 2x2 block for λ 1 = 2 \lambda_1 = 2 λ 1 = 2 and a 1x1 block for λ 2 = 1 \lambda_2 = 1 λ 2 = 1
A ˉ = J = [ 2 1 0 0 2 0 0 0 1 ] → P J = A P \bar A = J = \begin{bmatrix} 2 & 1 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 1 \end{bmatrix} \rightarrow \ P J = A P A ˉ = J = ⎣ ⎡ 2 0 0 1 2 0 0 0 1 ⎦ ⎤ → P J = A P
Now we will construct the chain of vectors, which will be used to construct the basis vectors.
x ∈ N ( Σ 1 2 ) but x ∉ N ( Σ 1 ) x \in N(\Sigma_1^2) \text{ but } x \notin N(\Sigma_1) x ∈ N ( Σ 1 2 ) but x ∈ / N ( Σ 1 )
Σ 1 2 x = 0 → [ 0 0 4 0 0 1 0 0 1 ] [ x 1 x 2 x 3 ] = [ 0 0 0 ] → x 3 = 0 \Sigma_1^2x = 0 \rightarrow \begin{bmatrix} 0 & 0 & 4 \\ 0 & 0 & 1 \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \\ x_3 \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ 0 \end{bmatrix} \rightarrow x_3 = 0 Σ 1 2 x = 0 → ⎣ ⎡ 0 0 0 0 0 0 4 1 1 ⎦ ⎤ ⎣ ⎡ x 1 x 2 x 3 ⎦ ⎤ = ⎣ ⎡ 0 0 0 ⎦ ⎤ → x 3 = 0
Σ 1 x ≠ 0 → [ − 1 1 − 1 − 1 1 2 0 0 − 1 ] [ x 1 x 2 0 ] ≠ [ 0 0 0 ] → x 1 ≠ x 2 \Sigma_1x \neq 0 \rightarrow \begin{bmatrix} -1 & 1 & -1 \\ -1 & 1 & 2 \\ 0 & 0 & -1 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \\ 0 \end{bmatrix} \neq \begin{bmatrix} 0 \\ 0 \\ 0 \end{bmatrix} \rightarrow x_1 \neq x_2 Σ 1 x = 0 → ⎣ ⎡ − 1 − 1 0 1 1 0 − 1 2 − 1 ⎦ ⎤ ⎣ ⎡ x 1 x 2 0 ⎦ ⎤ = ⎣ ⎡ 0 0 0 ⎦ ⎤ → x 1 = x 2
Select x = [ 1 0 0 ] x = \begin{bmatrix} 1 \\ 0 \\ 0 \end{bmatrix} x = ⎣ ⎡ 1 0 0 ⎦ ⎤ such that x ∈ N ( Σ 1 2 ) but x ∉ N ( Σ 1 ) x \in N(\Sigma_1^2) \text{ but } x \notin N(\Sigma_1) x ∈ N ( Σ 1 2 ) but x ∈ / N ( Σ 1 )
The chain of vectors for λ 1 = 2 \lambda_1 = 2 λ 1 = 2 is
{ Σ 1 x , x } = { [ − 1 − 1 0 ] , [ 1 0 0 ] } \Big \{ \ \Sigma_1x \ , \ x \ \Big \} = \Bigg \{ \begin{bmatrix} -1 \\ -1 \\ 0 \end{bmatrix}, \begin{bmatrix} 1 \\ 0 \\ 0 \end{bmatrix} \Bigg \} { Σ 1 x , x } = { ⎣ ⎡ − 1 − 1 0 ⎦ ⎤ , ⎣ ⎡ 1 0 0 ⎦ ⎤ }
┌──────────────────────────────────┬──────────┐
│ │ N(Σ₁²) │
│ - x └──────────┤
│ ┌───────────────┬──────────┤
│ │ │ N(Σ₁) │
│ │ └──────────┤
│ │ - Σ₁x │
└──────────────────┴──────────────────────────┘
Now we will construct the chain of vectors for λ 2 = 1 \lambda_2 = 1 λ 2 = 1
Σ 2 = A − I = [ 0 1 − 1 − 1 2 2 0 0 0 ] \Sigma_2 = A - I = \begin{bmatrix} 0 & 1 & -1 \\ -1 & 2 & 2 \\ 0 & 0 & 0 \end{bmatrix} Σ 2 = A − I = ⎣ ⎡ 0 − 1 0 1 2 0 − 1 2 0 ⎦ ⎤
dim N ( Σ 2 ) = 3 − 2 = 1 = r 2 \text{ dim }N(\Sigma_2) = 3-2 = 1 = r_2 dim N ( Σ 2 ) = 3 − 2 = 1 = r 2
y ∈ N ( Σ 2 ) y \in N(\Sigma_2) y ∈ N ( Σ 2 )
Σ 2 y = 0 → [ 0 1 − 1 − 1 2 2 0 0 0 ] [ y 1 y 2 y 3 ] = [ 0 0 0 ] → y = [ 4 1 1 ] \Sigma_2y = 0 \rightarrow \begin{bmatrix} 0 & 1 & -1 \\ -1 & 2 & 2 \\ 0 & 0 & 0 \end{bmatrix} \begin{bmatrix} y_1 \\ y_2 \\ y_3 \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \\ 0 \end{bmatrix} \rightarrow y = \begin{bmatrix} 4 \\ 1 \\ 1 \end{bmatrix} Σ 2 y = 0 → ⎣ ⎡ 0 − 1 0 1 2 0 − 1 2 0 ⎦ ⎤ ⎣ ⎡ y 1 y 2 y 3 ⎦ ⎤ = ⎣ ⎡ 0 0 0 ⎦ ⎤ → y = ⎣ ⎡ 4 1 1 ⎦ ⎤
The chain of vectors for λ 2 = 1 \lambda_2 = 1 λ 2 = 1 is
{ y } = { [ 4 1 1 ] } \Big \{y \Big \} = \Bigg \{\begin{bmatrix} 4 \\ 1 \\ 1 \end{bmatrix} \Bigg \} { y } = { ⎣ ⎡ 4 1 1 ⎦ ⎤ }
┌──────────────────────────────────┬──────────┐
│ │ N(Σ₂) │
│ └──────────┤
│ - y │
└─────────────────────────────────────────────┘
Now we will construct the change of basis matrix P P P , or B B B in the general case.
P = [ Σ 1 x ∣ x ∣ y ] = [ − 1 1 4 − 1 0 1 0 0 1 ] P = \begin{bmatrix} \ \ \Sigma_1x \ \ \ \ \vert & \ x & \vert \ \ \ \ \ y \ \ \end{bmatrix} = \begin{bmatrix} -1 & 1 & 4 \\ -1 & 0 & 1 \\ 0 & 0 & 1 \end{bmatrix} P = [ Σ 1 x ∣ x ∣ y ] = ⎣ ⎡ − 1 − 1 0 1 0 0 4 1 1 ⎦ ⎤
J = P − 1 A P = [ − 1 1 4 − 1 0 1 0 0 1 ] − 1 [ 1 1 − 1 − 1 3 2 0 0 1 ] [ − 1 1 4 − 1 0 1 0 0 1 ] = [ 2 1 0 0 2 0 0 0 1 ] J = P^{-1}AP \ = \ \begin{bmatrix} -1 & 1 & 4 \\ -1 & 0 & 1 \\ 0 & 0 & 1 \end{bmatrix}^{-1} \begin{bmatrix} 1 & 1 & -1 \\ -1 & 3 & 2 \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} -1 & 1 & 4 \\ -1 & 0 & 1 \\ 0 & 0 & 1 \end{bmatrix} = \begin{bmatrix} 2 & 1 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 1 \end{bmatrix} J = P − 1 A P = ⎣ ⎡ − 1 − 1 0 1 0 0 4 1 1 ⎦ ⎤ − 1 ⎣ ⎡ 1 − 1 0 1 3 0 − 1 2 1 ⎦ ⎤ ⎣ ⎡ − 1 − 1 0 1 0 0 4 1 1 ⎦ ⎤ = ⎣ ⎡ 2 0 0 1 2 0 0 0 1 ⎦ ⎤
A P = P A ˉ ⟹ [ 1 1 − 1 − 1 3 2 0 0 1 ] [ − 1 1 4 − 1 1 1 0 0 1 ] = [ − 1 1 4 − 1 1 1 0 0 1 ] [ 2 1 0 0 2 0 0 0 1 ] A P = P \bar A \implies \begin{bmatrix} 1 & 1 & -1 \\ -1 & 3 & 2 \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} -1 & 1 & 4 \\ -1 & 1 & 1 \\ 0 & 0 & 1 \end{bmatrix} = \begin{bmatrix} -1 & 1 & 4 \\ -1 & 1 & 1 \\ 0 & 0 & 1 \end{bmatrix} \begin{bmatrix} 2 & 1 & 0 \\ 0 & 2 & 0 \\ 0 & 0 & 1 \end{bmatrix} A P = P A ˉ ⟹ ⎣ ⎡ 1 − 1 0 1 3 0 − 1 2 1 ⎦ ⎤ ⎣ ⎡ − 1 − 1 0 1 1 0 4 1 1 ⎦ ⎤ = ⎣ ⎡ − 1 − 1 0 1 1 0 4 1 1 ⎦ ⎤ ⎣ ⎡ 2 0 0 1 2 0 0 0 1 ⎦ ⎤
Special Case I
A A A has a single eigenvalue λ i \lambda_i λ i and m i = r i m_i = r_i m i = r i
┌─────────────────────────┬──────┬───────┬──────┐
│ N(Σᵢ) - Σᵢᵐ⁻¹x │ │ ..... │ │
├─────────────────────────┘ │ │ │
│ N(Σᵢ²) - Σᵢᵐ⁻²x │ │ │
├────────────────────────────────┘.. │ │
│ . .. │ │
├────────────────────────────────────────┘ │
│ N(Σᵢᵐᶦ) - x │
└───────────────────────────────────────────────┘
For this case, we can see that m ( s ) = d ( s ) m(s) = d(s) m ( s ) = d ( s ) , which means that the Jordan canonical form will have a single block for λ \lambda λ , and the size of the block will be m i = r i m_i = r_i m i = r i . This is because the minimal polynomial is the product of the linear factors of the characteristic polynomial, and the characteristic polynomial is the product of the eigenvalues. Since there is only one eigenvalue, the minimal polynomial will be the same as the characteristic polynomial.
A ˉ = J = [ λ 1 0 ⋯ 0 0 λ 1 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 ⋯ λ ] \bar A = J = \begin{bmatrix} \lambda & 1 & 0 & \cdots & 0 \\ 0 & \lambda & 1 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & 0 & \cdots & \lambda \end{bmatrix} A ˉ = J = ⎣ ⎡ λ 0 ⋮ 0 1 λ ⋮ 0 0 1 ⋮ 0 ⋯ ⋯ ⋱ ⋯ 0 0 ⋮ λ ⎦ ⎤
consider the case: { Σ i m i − 1 x , Σ i m i − 2 x , ⋯ , Σ i x , x } \Big \{ \Sigma_i^{m_i-1}x, \Sigma_i^{m_i-2}x, \cdots, \Sigma_ix, x \Big \} { Σ i m i − 1 x , Σ i m i − 2 x , ⋯ , Σ i x , x } where it is sufficient to span N ( Σ i m i ) N(\Sigma_i^{m_i}) N ( Σ i m i ) , which is the same as N ( Σ i m i − 1 ) N(\Sigma_i^{m_i-1}) N ( Σ i m i − 1 ) , which is the same as N ( Σ i m i − 2 ) N(\Sigma_i^{m_i-2}) N ( Σ i m i − 2 ) , and so on. There are r i r_i r i vectors in this chain, which is the same as the dimension of N ( Σ i m i ) N(\Sigma_i^{m_i}) N ( Σ i m i ) .
Special Case II
A A A has a single eigenvalue λ i \lambda_i λ i and m i = 1 m_i = 1 m i = 1
┌───────────────────────────────────────────────┐
│ N(Σᵢ) - x₃ │
│ - x₂ │
│ - x₁ │
└───────────────────────────────────────────────┘
For this case, we can see that m ( s ) = ( s − λ i ) m(s) = (s-\lambda_i) m ( s ) = ( s − λ i ) , which means that the Jordan canonical form will have a single block for λ \lambda λ , and the size of the block will be m i = 1 m_i = 1 m i = 1 . This is because the minimal polynomial is the product of the linear factors of the characteristic polynomial, and the characteristic polynomial is the product of the eigenvalues. Since there is only one eigenvalue, the minimal polynomial will be the same as the characteristic polynomial. The number of chains will be equal to the dimension of N ( Σ i ) N(\Sigma_i) N ( Σ i ) , which is r i r_i r i .
dim N ( Σ i ) = r i \text{dim }N(\Sigma_i) = r_i dim N ( Σ i ) = r i
chain 1 = { x 1 } chain 2 = { x 2 } ⋮ chain r i = { x r i } \text{chain }1 = \Big \{ x_1 \Big \} \\
\text{chain }2 = \Big \{ x_2 \Big \} \\
\ \ \ \ \ \ \ \ \ \ \ \vdots \\
\text{chain }r_i = \Big \{ x_{r_i} \Big \} chain 1 = { x 1 } chain 2 = { x 2 } ⋮ chain r i = { x r i }
J = [ λ 0 0 ⋯ 0 0 λ 0 ⋯ 0 ⋮ ⋮ ⋮ ⋱ ⋮ 0 0 0 ⋯ λ ] J = \begin{bmatrix} \lambda & 0 & 0 & \cdots & 0 \\ 0 & \lambda & 0 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & 0 & \cdots & \lambda \end{bmatrix} J = ⎣ ⎡ λ 0 ⋮ 0 0 λ ⋮ 0 0 0 ⋮ 0 ⋯ ⋯ ⋱ ⋯ 0 0 ⋮ λ ⎦ ⎤
Example : Given m i = 6 m_i = 6 m i = 6 and r i = 8 r_i = 8 r i = 8 , Find the largest K value such that y ∈ N ( Σ i K ) y \in N(\Sigma_i^K) y ∈ N ( Σ i K ) but y ∉ N ( Σ i K − 1 ) y \notin N(\Sigma_i^{K-1}) y ∈ / N ( Σ i K − 1 )
┌─────────────────────────┬──────┬───────┬──────┬──────┬───────┐
│ N(Σᵢ) -Σᵢy - Σᵢ⁵ │ │ │ │ │ │
├─────────────────────────┘ │ │ │ │ │
│ N(Σᵢ²) -y - Σᵢ⁴x │ │ │ │ │
├────────────────────────────────┘ │ │ │ │
│ N(Σᵢ³) - Σᵢ³x │ │ │ │
├────────────────────────────────────────┘ │ │ │
│ N(Σᵢ⁴) - Σᵢ²x │ │ │
├───────────────────────────────────────────────┘ │ │
│ N(Σᵢ⁵) - Σᵢx │ │
├──────────────────────────────────────────────────────┘ │
│ N(Σᵢ⁶) - x │
└──────────────────────────────────────────────────────────────┘
First chain { Σ i 5 x , Σ i 4 x , Σ i 3 x , Σ i 2 x , Σ i x , x } \Big \{ \Sigma_i^5x, \Sigma_i^4x, \Sigma_i^3x, \Sigma_i^2x, \Sigma_ix, x \Big \} { Σ i 5 x , Σ i 4 x , Σ i 3 x , Σ i 2 x , Σ i x , x } has r i r_i r i vectors, which is the same as the dimension of N ( Σ i 6 ) N(\Sigma_i^6) N ( Σ i 6 )
Second chain { Σ i y , y } \Big \{ \Sigma_i y, y \Big \} { Σ i y , y } has 2 vectors, which is the same as the dimension of N ( Σ i 1 ) N(\Sigma_i^1) N ( Σ i 1 )
B = [ Σ i 5 x Σ i 4 x Σ i 3 x Σ i 2 x Σ i x x Σ i y y ] B = \begin{bmatrix} \Sigma_i^5x & \Sigma_i^4x & \Sigma_i^3x & \Sigma_i^2x & \Sigma_ix & x & \Sigma_iy & y \end{bmatrix} B = [ Σ i 5 x Σ i 4 x Σ i 3 x Σ i 2 x Σ i x x Σ i y y ]
J = [ λ 1 0 0 0 0 0 0 0 λ 1 0 0 0 0 0 0 0 λ 1 0 0 0 0 0 0 0 λ 1 0 0 0 0 0 0 0 λ 1 0 0 0 0 0 0 0 λ 0 0 0 0 0 0 0 0 λ 1 0 0 0 0 0 0 0 λ ] J = \begin{bmatrix} \lambda & 1 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & \lambda & 1 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & \lambda & 1 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & \lambda & 1 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & \lambda & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & \lambda & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 & \lambda & 1 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 & \lambda \end{bmatrix} J = ⎣ ⎡ λ 0 0 0 0 0 0 0 1 λ 0 0 0 0 0 0 0 1 λ 0 0 0 0 0 0 0 1 λ 0 0 0 0 0 0 0 1 λ 0 0 0 0 0 0 0 1 λ 0 0 0 0 0 0 0 0 λ 0 0 0 0 0 0 0 1 λ ⎦ ⎤
Example : Let A ∈ C 6 × 6 A \in \mathbb{C}^{6\times 6} A ∈ C 6 × 6 with a single eigenvalue. m ( s ) = ( s − λ i ) 3 m(s)=(s-\lambda_i)^3 m ( s ) = ( s − λ i ) 3 , dim N ( Σ i ) = 2 \text{dim }N(\Sigma_i) = 2 dim N ( Σ i ) = 2 , Find the Jordan canonical form of A A A and construct the change of basis matrix B B B .
Solution :
d ( s ) = ( s − λ i ) 6 → dim ( v ) = 6 single eigenvalue λ i d(s) = (s - \lambda_i )^6 \rightarrow \text{dim }(v) = 6 \text{ single eigenvalue } \lambda_i d ( s ) = ( s − λ i ) 6 → dim ( v ) = 6 single eigenvalue λ i
dim N ( Σ i ) = 2 ⟹ There are 2 Jordan blocks or chains \text{dim }N(\Sigma_i) = 2 \implies \text{ There are 2 Jordan blocks or chains } dim N ( Σ i ) = 2 ⟹ There are 2 Jordan blocks or chains
m ( s ) = ( s − λ i ) 3 ⟹ The size of largest Jordan block is 3 × 3 m(s) = (s-\lambda_i)^3 \implies \text{ The size of largest Jordan block is } 3 \times 3 m ( s ) = ( s − λ i ) 3 ⟹ The size of largest Jordan block is 3 × 3
Then we would have 2 Jordan blocks of size 3x3. The change of basis matrix B B B would be composed of the basis vectors.
A ˉ = J = [ λ 1 0 0 0 0 0 λ 1 0 0 0 0 0 λ 0 0 0 0 0 0 λ 1 0 0 0 0 0 λ 1 0 0 0 0 0 λ ] \bar A = J = \begin{bmatrix} \lambda & 1 & 0 & 0 & 0 & 0 \\ 0 & \lambda & 1 & 0 & 0 & 0 \\ 0 & 0 & \lambda & 0 & 0 & 0 \\ 0 & 0 & 0 & \lambda & 1 & 0 \\ 0 & 0 & 0 & 0 & \lambda & 1 \\ 0 & 0 & 0 & 0 & 0 & \lambda \end{bmatrix} A ˉ = J = ⎣ ⎡ λ 0 0 0 0 0 1 λ 0 0 0 0 0 1 λ 0 0 0 0 0 0 λ 0 0 0 0 0 1 λ 0 0 0 0 0 1 λ ⎦ ⎤
Jordan chains for λ i \lambda_i λ i are { Σ i 2 x , Σ i x , x } \Big \{ \Sigma_i^2x, \Sigma_ix, x \Big \} { Σ i 2 x , Σ i x , x } and { Σ i 2 y , Σ i y , y } \Big \{ \Sigma_i^2y, \Sigma_iy, y \Big \} { Σ i 2 y , Σ i y , y }
B = [ Σ i 2 x Σ i x x Σ i 2 y Σ i y y ] B = \begin{bmatrix} \Sigma_i^2x & \Sigma_ix & x & \Sigma_i^2y & \Sigma_iy & y \end{bmatrix} B = [ Σ i 2 x Σ i x x Σ i 2 y Σ i y y ]
Try to find x x x and y y y such that x ∈ N ( Σ i 2 ) x \in N(\Sigma_i^2) x ∈ N ( Σ i 2 ) but x ∉ N ( Σ i ) x \notin N(\Sigma_i) x ∈ / N ( Σ i ) and y ∈ N ( Σ i ) y \in N(\Sigma_i) y ∈ N ( Σ i )
ASK FOR HELP FOR NEXT
Example : Let A = [ 1 1 0 0 0 1 0 0 0 0 1 0 − 1 1 0 2 ] A = \begin{bmatrix} 1 & 1 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ -1 & 1 & 0 & 2 \end{bmatrix} A = ⎣ ⎡ 1 0 0 − 1 1 1 0 1 0 0 1 0 0 0 0 2 ⎦ ⎤ Find the characteristic polynomial, minimal polynomial, Jordan canonical form and the basis for A A A .
Solution :
d ( s ) = det ( s I − A ) = ( s − 1 ) 3 ( s − 2 ) d(s) = \text{det}(sI - A) = (s-1)^3(s-2) d ( s ) = det ( s I − A ) = ( s − 1 ) 3 ( s − 2 )
For the minimal polynomial, we need to check the null space of A − λ i I A - \lambda_iI A − λ i I for each eigenvalue.
Σ 1 = A − I = [ 0 1 0 0 0 0 0 0 0 0 0 0 − 1 1 0 1 ] dim N ( Σ 1 ) = 4 − 2 = 2 ≠ r 1 \Sigma_1 = \ \ A - I \ \ = \begin{bmatrix} 0 & 1 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ -1 & 1 & 0 & 1 \end{bmatrix}
\text{dim }N(\Sigma_1) = 4-2 = 2 \neq r_1 Σ 1 = A − I = ⎣ ⎡ 0 0 0 − 1 1 0 0 1 0 0 0 0 0 0 0 1 ⎦ ⎤ dim N ( Σ 1 ) = 4 − 2 = 2 = r 1
Σ 1 2 = ( A − I ) 2 = [ 0 0 0 0 0 0 0 0 0 0 0 0 − 1 0 0 1 ] dim N ( Σ 1 2 ) = 4 − 1 = 3 = r 1 \Sigma_1^2 = (A - I)^2 = \begin{bmatrix} 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 \\ -1 & 0 & 0 & 1 \end{bmatrix}
\text{dim }N(\Sigma_1^2) = 4-1 = 3 = r_1 Σ 1 2 = ( A − I ) 2 = ⎣ ⎡ 0 0 0 − 1 0 0 0 0 0 0 0 0 0 0 0 1 ⎦ ⎤ dim N ( Σ 1 2 ) = 4 − 1 = 3 = r 1
Σ 2 = A − 2 I = [ − 1 1 0 0 0 − 1 0 0 0 0 − 1 0 − 1 1 0 0 ] dim N ( Σ 2 ) = 4 − 1 = 3 = r 2 \Sigma_2 = \ \ A - 2I \ = \begin{bmatrix} -1 & 1 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & -1 & 0 \\ -1 & 1 & 0 & 0 \end{bmatrix}
\text{dim }N(\Sigma_2) = 4-1 = 3 = r_2 Σ 2 = A − 2 I = ⎣ ⎡ − 1 0 0 − 1 1 − 1 0 1 0 0 − 1 0 0 0 0 0 ⎦ ⎤ dim N ( Σ 2 ) = 4 − 1 = 3 = r 2
m ( s ) = ( s − 1 ) 2 ( s − 2 ) m(s) = (s-1)^2(s-2) m ( s ) = ( s − 1 ) 2 ( s − 2 ) Largest Jordan block is 2x2
A ˉ = J = [ 1 1 0 0 0 1 0 0 0 0 1 0 0 0 0 2 ] \bar A = J = \begin{bmatrix} 1 & 1 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 2 \end{bmatrix} A ˉ = J = ⎣ ⎡ 1 0 0 0 1 1 0 0 0 0 1 0 0 0 0 2 ⎦ ⎤
Jordan chains for λ 1 = 1 \lambda_1 = 1 λ 1 = 1 are { Σ 1 x , x } \Big \{ \Sigma_1x, x \Big \} { Σ 1 x , x }
choose x = [ 1 1 0 1 ] x = \begin{bmatrix} 1 \\ 1 \\ 0 \\ 1 \end{bmatrix} x = ⎣ ⎡ 1 1 0 1 ⎦ ⎤
Then { Σ 1 x , x } = { [ 1 0 0 1 ] , [ 1 1 0 1 ] } \Big \{ \Sigma_1x, x \Big \} = \Big \{ \begin{bmatrix} 1 \\ 0 \\ 0 \\ 1 \end{bmatrix}, \begin{bmatrix} 1 \\ 1 \\ 0 \\ 1 \end{bmatrix} \Big \} { Σ 1 x , x } = { ⎣ ⎡ 1 0 0 1 ⎦ ⎤ , ⎣ ⎡ 1 1 0 1 ⎦ ⎤ }
choose y = [ 0 0 1 0 ] y = \begin{bmatrix} 0 \\ 0 \\ 1 \\ 0 \end{bmatrix} y = ⎣ ⎡ 0 0 1 0 ⎦ ⎤
Then { y } = { [ 0 0 1 0 ] } \Big \{y \Big \} = \Big \{ \begin{bmatrix} 0 \\ 0 \\ 1 \\ 0 \end{bmatrix} \Big \} { y } = { ⎣ ⎡ 0 0 1 0 ⎦ ⎤ }
choose z = [ 0 0 0 1 ] z = \begin{bmatrix} 0 \\ 0 \\ 0 \\ 1 \end{bmatrix} z = ⎣ ⎡ 0 0 0 1 ⎦ ⎤
Then { z } = { [ 0 0 0 1 ] } \Big \{z \Big \} = \Big \{ \begin{bmatrix} 0 \\ 0 \\ 0 \\ 1 \end{bmatrix} \Big \} { z } = { ⎣ ⎡ 0 0 0 1 ⎦ ⎤ }
B = [ Σ 1 x x y z ] = [ 1 1 0 0 0 1 0 0 0 0 1 0 1 1 0 1 ] B = \begin{bmatrix} \Sigma_1x & x & y & z \end{bmatrix} = \begin{bmatrix} 1 & 1 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 1 & 1 & 0 & 1 \end{bmatrix} B = [ Σ 1 x x y z ] = ⎣ ⎡ 1 0 0 1 1 1 0 1 0 0 1 0 0 0 0 1 ⎦ ⎤
#EE501 - Linear Systems Theory at METU